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This  papar  presents  a  pr  alia  in  ary  study  on 
using  Mathanatical  Morphology  to  raprasant  and 
coda  a  binary  or  a  grey-tone  image  by  parts  of 
its  skalaton*  a  thinnad  varsion  of  tba  image.  An 
image  can  ba  uniquely  decomposed  into  skeleton 
components/  and  then  reconstructed  by  dilating 
these  components.  Since*  for  a  certain  category 
of  imagery*  the  skeleton  components  possess  a 
lower  entropy  than  the  original  image*  a  run- 
length  or  entropy  coding  scheme  can  be  used  to 
achieve  representation  or  transmission  of  tba 
image  at  a  lower  information  rate  than  originally 
required. 

Mathematical  Morphology*  as  a  method  for 
image  analysis*  was  introduced  by  Matheron  and 
Berra  P  ]  •  Its  purpose  is  the  quantitative  des¬ 
cription  of  geometrical  structures.  TO  extract 
information  from  an  image  object*  Morphology 
•hits*  it  first  with  a  "structuring  element." 
The  Interaction  with  the  structuring  element 
transforms  the  object  and  reduces  it  to  a  sort  of 
caricature  which  is  more  expressive  than  the 
actual  initial  phenomenon. 

The  most  fundmsental  morphological  transfor¬ 
mations  are  erosion  and  dilationi  Let  X  denote  a 
set  in  the  continuous  or  digital  2-D  Euclidean 
space  representing  a  binary  analog  or  digital 
image  object.  Then  X°  {complement  of  X)  denotes 
the  image  background.  Let  B  be  the  structuring 
element*  which  is  another  set  with  a  simple  geo¬ 
metrical  shape,  and  denote  by  Bx  the  translate  of 
B  whose  center  is  situated  at  the  point  x. 
Erosion  of  X  by  B  is  the  set  of  all  points  x  such 
that  B^  is  included  in  X  (see  Fig.  1).  Symboli¬ 
cally, 

X0B-{xt  »aCX)  (1) 

Dilation  of  X  by  B  is  the  set  of  all  points  x 
such  that  Ss  "hits"  Xy  i.e.  has  a  non1  empty 
intersection  with  X.  Symbolically, 

(2> 

Fig.  1  shows  the  erosion  end  dilation  of  e  set  X 
by  a  disk  I.  This  figure  illustrates  that  ire* 


sion  is  a  shrinking  operation  and  dilation  is  en 
expanding  operation.  Erosion  end  dilation  are 
dual  operations  w.r •  to  complementations  Eroding 
X  is  equivalent  to  taking  the  complement  of  the 
dilation  of  X°.  If  we  erode  X  by  B  end  then 
dilate  the  set  X0B  by  B,  we  do  not  reoover  X. 
Me  reconstitute  only  e  part  of  X  which  la  simpler 
end  has  lees  details.  It  may  be  considered  as 
that  part  which  is  most  essential  morphological¬ 
ly.  Me  call  this  new  set  the  opening  of  X  w.r. 
to  Bi 

*«’(X0B)  ©b  (3) 

The  opening  is  the  domain  swept  out  by  all  the 
translates  of  r*  which  are  included  in  X.  This 
operation  amooths  the  contours  of  X,  outs  the 
narrow  Isthmuses,  suppresses  the  small  islands 
and  the  sharp  capes  of  X. 


Although  the  above  operations  appear  super¬ 
ficially  simple,  we  can  perform  an  enormous 
variety  of  image  processing  and  image  understand¬ 
ing  tasks  just  by  combining  erosions  sad 
dilations,  as  is  well  developed  in  (1]. 


The  skeleton  is  a  topologically  equivalent 
thinned  version  of  the  image.  It  can  be  obtained 
from  morphological  transformations  which  empha¬ 
sise  features  of  the  object  associated  with  its 
connectivity.  In  tbs  2-D  continuous  speot  it  is 
defined  as  follows  t  Let  rB^  denote  the  disk  of 
radius  r  cantered  at  the  point  x.  Let  sc  (X) 
denote  the  set  of  the  oeaters  of  the  disks  rB_ 
such  thatt  i)  rB  is  the  maxima  disk  oentareo 
at  x  md  contained  in  the  object  X,  and  li)  tba 
disk  rB  intarseets  the  boundary  of  X  at  two  or 
more  different  places.  Then,  the  skeleton  B(X) 
of  x  is  defined  as  the  set  of  the  oeaters  of  the 
max  in  us  disks  laser  ihable  la  X,  sad  is  s  carica¬ 
ture  containing  information  about  the  shape,  site 
and  orientation  of  X.  Borne  examples  of  skelofeoos 
are  shown  in  Fig.  2.  The  skeleton  B  (X)  eaa  be 
obtained  from  the  set  onion  of  s.  (X)  (Lantuejoul 
(21) » 
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where  *U*  ("/")  represents  Mt  union  (dif¬ 

ference)  #  and  dr  in  the  infinitesimal  Mall 
radius. 


Although  tha  skalaton  la  not  a  well-dig ite^ 
llsable  notion#  Serra  (1)  gives  an  algorithm  for 
tha  skalaton  of  digital  binary  images  sampled  on 
a  hax agonal  grid. 


Our  rasaarch  was  focused  on  thraa  arsast 
obtaining  algorithms  for  skalatonising  digital 
binary  Images  on  a  ractangular  grid;  using  parts 
of  tha  skalaton  to  coda  tba  image;  and  extending 
tha  above  ideas  to  gray- tone  images. 


Let  R  denote  tha  unit-sisa  square  of  a  rec¬ 
tangular  grid  (sea  Fig.  3)  which  is  a  square  of 
3x3  pixels#  and  let  nR  denote  tha  square  R  Mini¬ 
fied  n  tines  which  gives  s  square  of  (2n+1)  x 
(2n+1)  pixels.  Than  a  digital  algorithm  for  8{X) 
of  a  rectangularly  swpled  image  object  X  is 


S(X)  -  0  .#(X)  -  D  [x©n*)/(X©n*)R]  (5) 


n-0 


n«0 


lq.  (5)  says  that  tha  skalaton  subsets  s  (X)  torn 
a  partition  of  S(X).  Thus#  8(X)  is  obtained  by 
successively  eroding  X  by  nR#  and  than  keeping 
from  every  eroded  set  (X  (£)nX)  those  parts  only 
which  oonsist  of  angular  points  and  lines  without 
thickness;  these  parts  are  tha  only  ones  remain¬ 
ing  after  tha  set  difference  between  <X0nR)  snd 
its  opening  (X0nR)Jt*  Tha  maxima  site 
indicates  tha  square  of  maximal  sisa  after  which 
a  further  erosion  erodes  X  down  to  the  empty  set. 


Mow#  the  image  X  can  be  exactly  reconstruct¬ 
ed  by  dilating  the  subsets  of  its  skeleton  by 
squares  of  corresponding  sisa  and  taking  their 
unions 


0  [•„<!>  ©«*J 

n*0 


<«) 


■qs.  (5)  and  (ft)  imply  that  the  datum  of  the 
Image  (set  X)  is  equivalent  to  that  of  its  skele¬ 
ton  ft  (X)  together  with  the  else  "n*  of  the 
maximum  square  associated  with  each  point  of 
8(X) •  In  rig.  4#  proceeding  from  left  to  right 
columns#  we  show  an  example  of  an  image  object  X 
and  its  erosions  (X  Qnl)  #  the  openings  of  these 


erosions#  the  skeleton  subsets  s_(X)#  the  dilated 


the  composition  of  the  skeleton  ft(X)  as 
the  union  of  the  skeleton  subsets#  and  finally 
the  reconstruction  of  X  as  the  union  of  the  di¬ 
lated  skeleton  subsets. 


According  to  Shannon's  theory  of  dieerete 
coding  (31  we  oonslder  the  digitised 
images  as  sample  functions  of  a  3-D  stochastic 
process  characterised  by  joint  probability  dis¬ 
tributions  of  all  orders.  In  practice  we  measure 
histograms  instead  of  probability  distribu¬ 
tions.  Consider  a  1-0  or  2-D  block  of 


consecutive  pixels  x1#x2#»..#x||#  where  x#  can  be 
either  1  or  0  according  to  whether  x*  belongs  to 
the  image  object  X  or  its  background  x9  respec¬ 
tively.  Let  p(x1#x2#...#xv)  be  the  Mtb-order 
joint  probability  of  these  M  pixels.  Then  the 
Nth-order  joint  entropy  (in  bits/pixel)  of  the 
binary  image  XDX*  is  defined  as 


X^  #  •  •  •  #Xjj 


•  log2p(x1,...#x||)  (7) 


As  is  well  known#  ^  is  a  nonincreasing  function 
of  N  and  the  limit  as  M  ♦  •  is  the  entropy of  the 
stochastic  source.  If  we  oonslder  the  2*  dif¬ 
ferent  blocks  of  M  pixels  each  as  our  messages# 
we  can  employ  Huffman  coding  or  other  subopt  la  urn 
coding  procedures  [4]  to  achieve  transmission 
rates  very  cloee  to  these  Nth-order  entropies. 
Thus#  hereafter  we  will  be  referring  to  these 
Nth-order  joint  entropies  of  binary  Images  as 
their  achievable  transmission  rates. 

Since  every  skeleton  subset  sR(X)  is  •  much 
thinner  binary  image  than  X#  then  its  Mtb-order 
entropy#  denoted  by  H,|(sR)#  will  be  much  lower 
than  Bp  (A).  And  there  might  be  cases  where 

^max 

l  ■-(.  )  «  *,<*>  (•) 

n«0  n  " 

Thus#  to  transmit  s  (X)  we  need  approximately 
hm(sr)  bits/pixel.  In  addition  to  the  sta  of  all 
Bg (s_)  we  need  information  about  which 

can  be  taken  into  account  with  the  trivial  amount 
of  log2 (M/2)  bits#  for  a  binary  image  of  nm 
pixels. 

When  (t)  holds#  we  can  transmit  all  the 
skeleton  subsets  of  X  Independently  at  a  total 
rate  less  than  the  entropy  of  the  original  image# 
and  fully  reconstruct  X  without  error  as  Bq.  (ft) 
indicates. 

A  further  reduction  in  information  rate  can 
be  achieved  by  using  not  all  but  only  same  of  the 
skeleton  subsets  to  reconstruct  openings  (smooth¬ 
ed  versions)  of  the  original  images 


• 0  ivx)  ©■*)  <•> 

n*k 

That  is#  if  in  the  imion  of  the  skeleton  subsets 
we  emit  the  first  k  subsets  (n*0#  1# .  •  #k-1 )  #  we 
reconstruct  the  opening  of  X  w.r.  to  kR.  The 
larger  the  k#  the  fewer  subsets  we  transmit#  the 
more  we  reduce  the  information  rate#  but  the 
smoother  is  the  version  X  that  we  reconstruct. 
As  shown  in  the  example  of  Pig  4#  for  N*4#  the 
original  image  X  has  an  entropy  of  0.34  bits/pix¬ 
el.  Xf  we  use  all  the  skeleton  subsets  we 


reconstruct  X  perfectly  st  s  rate  of  approxi¬ 
mately  0,16  bits/pixel.  If  we  desire  to 
reconstruct  only  the  openings  Xg  or  X^,  we  omit 
the  first  one  or  two  skeleton  subsets  and  thus  we 
need  approximately  0*16  or  0.14  bits/pixel  re¬ 
spectively.  Table  1  illustrates  that  more  infor¬ 
matively. 


BSU  1 

Mtb-Order  Atroples  (bits/pimel)  of  a  Adatom 
Baooostr noted  image  sod  Xta  Openings. 


M 

Image 

1 

2 

4 

8 

X 

0.47 

0.22 

0.18 

0.15 

** 

0.22 

0.19 

0.16 

0.13 

*2R 

0.20 

0.17 

0.14 

0.10 

0.07 

0.06 

0.05 

0.03 

The  first-. 

second-. 

fourth-  and 

eighth-order 

entropies  of 

the  original  binary 

image 

without 

skeleton  encoding  are  0.79c  0.50  ,  0.34  and  0.23 
respectively.  Thus,  as  shown  in  Table  1,  the  sum 
of  the  entropies  of  all  or  some  of  the  skeleton 
subsets  is  smaller  than  the  entropies  of  the 
original  unencoded  image. 


In  grey-tone  Morphology  111  the  binary  ero¬ 
sions  and  dilations  are  replaced  by  "min*  and 
"max”  operators  respectively.  Consider  s  nonne¬ 
gative  bounded  function  f(i,j)  representing  the 
intensity  of  a  sampled  grey-tone  image.  Let 
0<ttl,))<m  for  every  integer  pair  (i,j)  in  the 
image  support.  All  the  sero-valued  image  samples 
will  belong  by  convention  to  the  background  of 
the  image  object.  Erosion  or  dilation  of  the 
function  f  by  the  2-D  structuring  element  Ik  is 
defined  [1]  as 


[f  ©*](i#j)-min[f (r,s)«  ]  (10a) 

[*©r](l,j)-M*[f(r,.)«  (r,.)dt(1  .j)  (10b) 


of  (11)  will  be  aonnegative  functions.  Similarly 
as  in  Bq.  (6)  or  (9) ,  the  function  f  or  its 
openings  can  be  reconstructed  by  summing  alge¬ 
braically  all  or  some  of  the  skeleton 
subfunctions  sB(f)  dilated  by  Ol. 

The  implications  and  the  coding  efficiency 
of  the  skeleton  of  the  Image  function  f  in  terms 
of  entropy  considerations  are  still  under  inves¬ 
tigation. 


The  results  of  this  study  indicate  that  a 
digital  binary  image  can  be  uniquely  decomposed 
into  its  skeleton  and  the  maximum  inner ibable 
squares,  and  uniquely  reconstructed  from  its 
skeleton.  The  skeleton  provides  useful  informa¬ 
tion  about  the  shape,  else  and  orientation  of  an 
image.  Por  certain  categories  of  images  the 
total  entropy  of  the  skeleton  subsets  Is  lower 
than  the  entropy  of  the  original  images.  Origi¬ 
nal  1  bit/pixel  test  images  of  Irregularly  and 
regularly  shaped  objects  were  reconstructed  with¬ 
out  error  by  their  full  skeleton  st  information 
rates  of  *0.20  bits/pixel.  Smoothed  versions  of 
these  images  required  rates  of  only  ■ti.lS 
bits/pixel.  Finally,  by  using  min/max  operations 
instead  of  binary  eros ions/d ilat ions,  these  ideas 
can  be  extended  to  grey-tone  images. 
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where  R. ^  denotes  the  square  R  centered  at  the 
pixel  (l,)l.  The  opening  L  of  f  w.r.  to  R  is 
defined  as  an  erosion  followA  by  dilation. 

Me  provide  now  a  digital  algorithm  for  the 
skeleton  0(f)  of  f  which  will  be  the  nonnegative 
functions 

nrrr  *Vax 

■  (f)  -  l  .(f)  -  l  [ (f ©n*)- (f ©“*). ] 

n-0  B  a*>  *  (11) 

Bq.  (11)  is  s  direct  transposition  of  Bq.  (5) 
where  we  replaced  the  binary  set  union/diff erenow 
by  an  algebraic  addition/subtraction.  Because 
the  opening  is  an  anti-extensive  operation 
<fjl<f),  the  skeleton  subf unctions  in  the  brackets 


EROSION 


DILATION 


o 

X  ©  B  B 

Dilation  a  tit  X  by  B  ( after  [  2  ]  }. 


X©B 


Figure  2  -  Examples  of  Skeletons  (after  l  1  3). 


Figure  3-  The  3x3  pixels  square  R 
on  a  rectangular  grid. 
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Figure  4  -  Step  by  step  decomposition  and 
reconstruction  of  an  image  object  X  by  the 
components  of  its  skeleton  S(X)  : 

(a)  size  —  n  -  of  the  structuring  square  nR 

(b)  eroded  sets  (  X@nR  ) 

to  openings  of  the  eroded  sets  ( X©nR  )R 
(d)  skeleton  subsets  sn(X) 

<e)  dilated  skeleton  subsets  sn(X)®nR 

(f)  set  union  of  skeleton  subsets  $k(X)  for 

k-7,6 . n-Vn 

(g)  set  union  of  dilated  subsets  sk(X)®kR 
for  k  ■  7, 6, ,  n*1,  n  ,  which  glees  the 
opting  X„R 

(h)  twin  of  the  witropiet  H.( «.  )  of  the 
•Bwi  *k^X) ,  k  ■  7, ....  n ,  which  ere 
'•Wind  to  rcoonttruct  Ac  opening  X  n 
of  the  original  object  X 


